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Abstract 



A recently proposed approach for performing electronic-structure calculations 
on crystalline insulators in terms of localized orthogonal orbitals is applied to 
the oxides of lithium and sodium, L^O and Na20. Cohesive energies, lattice 
constants and bulk moduli of the aforementioned systems are determined 
at the Hartree-Fock level, and the corresponding values are shown to be in 
excellent agreement with the values obtained by a traditional Bloch-orbital- 
based Hartree-Fock approach. The present Wannier-function-based approach 
is expected to be advantageous in the treatment of electron-correlation effects 
in an infinite solid by conventional quantum-chemical methods. 
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I. INTRODUCTION 



A typical quantum-chemical investigation of a system, employing a wave-function-based 
approach, begins with a Hartree-Fock (HF) calculation which provides an initial mean-field 
description of the system. Then, if the need arises, it is improved systematically — the 
influence of electron correlations is included by considering virtual excitations from the HF 
wave function. However, in such a scheme the computational effort has a rather unfavor- 
able dependence on the size of the system and in the best-case scenario it scales roughly 
as iV 5 , where N is the number of atoms in the system. (Only recently, attempts towards 
achieving linear scaling at the HF levefi and towards exploiting the local nature of correla- 
tion effects are being madeH.) Because of this proliferation in the computational effort, the 
field of quantum chemistry is generally regarded as a science of finite systems with most of 
its applications limited to small or medium-sized molecules. Therefore, a naive extension 
of this quantum-chemical scheme to treat a system such as a three-dimensional crystalline 
solid is bound to run into problems. On the one hand a crystalline solid, for all practical 
purposes, is an infinite system. It essentially has an infinite number of electrons as well 
as an infinite spatial extension compared to the dimensions of a unit cell. On the other 
hand, the translational symmetry force its orbitals to have a delocalized itinerant charac- 
ter. Clearly, with such crystal orbitals, known as Bloch orbitals, it can be problematic to 
sum the contributions of the electron repulsion part of the Hamiltonian to convergence. In 
Bloch-orbital-based theories, the problems of infinite size and delocalized orbitals have been 
overcome at the Hartree-Fock level, nowadays, by adopting suitable procedures of integration 
over the Brillouin zone and lattice sums, as is done, for example, in the program CRYS- 
TAL!. However, inclusion of electron-correlation effects by going beyond the Hartree-Fock 
level within a traditional quantum-chemical scheme (such as the configuration-interaction or 
coupled-cluster approaches) appears to be fraught with problems if one uses Bloch orbitals 
as a basis set. Complications will arise both due to the itinerant nature of the Bloch orbitals, 
as well as the infinite number of virtual states spanning the Brillouin zone that one will have 
to take into account, in order to compute quantities such as the correlation energy. 

If one abandons the built-in translational symmetry associated with the Bloch orbitals, 
one can alternatively describe the electrons in a solid as being localized entities associated 
with the atoms or bonds constituting it. Such a localized description of electrons in a solid 
is much easier to understand for a chemist who is usually thinking in terms of bonds. Unlike 
the Bloch-orbital-based approach, a localized-orbital-oriented description of a system is not 
unambiguous, e.g. in the sense that one can choose either an orthogonal or a nonorthogonal 
set of localized orbitals. Nonorthogonal orbitals, in general, are better localized and, in 
principle, are also capable of describing the metallic systems where the electrons are indeed 
quite itinerant. Recently Bellaiche and Leyy@ have presented an approach, along with its 
application to crystalline LiH, which allows the determination of nonorthogonal localized 
Hartree-Fock orbitals of an infinite solid, from a set of finite-cluster calculations. However, 
with nonorthogonal orbitals a conventional quantum-chemical treatment of electron corre- 
lation such as the one based on the configuration-interaction approach is considerably more 
complicated. On the other hand, the orthogonal orbitals, although less localized than their 
nonorthogonal counterparts, offer the possibility of straightforward inclusion of correlation 
effects. Such an alternative approach to the description of crystalline electrons in terms 
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of localized orthogonal orbitals was proposed by Wannier long ag oB, and therefore, such 
orbitals in the condensed-matter community have come to be known as Wannier functions. 
It should be intuitively obvious that it will be very difficult to describe metallic electrons 
in terms of Wannier functions because the orthogonality requirement will force the orbitals 
to extend far away from their centers, leading to their eventual derealization!. However, 
for systems with a band gap such as semiconductors and insulators it should be possible, 
in principle, to satisfy both the orthogonality as well as the localization requirements as 
embodied in Wannier functions!. 

Therefore, it should be clear from the preceding discussion, that if one wants to study the 
influence of electron correlations on an infinite solid using a traditional wave-function-based 
approach, it is mandatory that one first obtains a HF representation of the system in terms 
of Wannier functions. This goal can be achieved by first performing a HF calculation for the 
infinite crystal using a Bloch-orbital-based approach, and then localizing the Bloch orbitals 
using one of the many available localization schemes! to obtain the Wannier functions!. 
However, in our previous papers§0 we had presented an alternative scheme which allows 
a direct determination of Hartree-Fock Wannier functions of a crystalline insulator. In the 
first paper (henceforth to be referred to as I), we had presented an outline of the formalism 
and used it to compute the total energy per unit cell of crystalline LiH within an LCAO 
scheme. In the second paper (to be called II), we gave a rigorous derivation of the approach 
by minimizing the total energy of an infinite crystal, and demonstrated its theoretical equiv- 
alence to a traditional Bloch-orbital-based HF approach. Additionally in II we applied the 
approach to obtain the HF Wannier functions of LiF and LiCl crystals which were subse- 
quently used to compute the X-ray structure factors and directional Compton profiles. In all 
the cases it was demonstrated that our results were in excellent agreement with equivalent 
calculations performed using the CRYSTAL program, which employs a Bloch-orbital-based 
HF approach. All the details pertaining to the computer implementation of our LCAO-HF 
program were also presented in 110. Recently, we extended this approach to study the band 
structure of insulating crystals and applied it to the case of NaCM We feel, however, that 
it is important for us to study the efficacy of the present approach in a variety of insulat- 
ing systems. With that goal in mind, in the present paper, we present the results of our 
Wannier-function-based HF calculations on crystalline Li 2 and Na20. Both these crystals 
are ionic in nature with the valence Wannier functions residing on O . However, O ion 
is not stable in its free form and becomes stabilized in the solid state only due to the crystal 
field. Therefore, we believe that obtaining the valence Wannier functions of compounds in- 
volving O is a nontrivial application of our approach. Moreover, this application involves 
rather diffuse basis functions requiring a significant improvement of the accuracy of lattice 
sums with respect to our calculations described in II. 

The remainder of our paper is organized as follows. In section |T[ we briefly review the 
theoretical formalism. Then in section |T| we present the results of our calculations for the 
Li 2 and Na 2 crystals. These results include cohesive energies, lattice constants and bulk 
moduli of the two compounds at the Hartree-Fock level. Our results are also shown to be in 
excellent agreement with those obtained using the CRYSTAL program. Finally, in section 
[TV| we present our conclusions. 
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II. THEORY 



In this section we briefly review the theory behind our approach. For more detailed 
derivations we refer the reader to papers I and IlSB. The essential idea behind the approach 
is that in order to describe a perfect crystalline insulator in terms of Wannier functions, one 
only needs to obtain the Wannier-type orbitals localized in any one of its unit cells — the 
orbitals localized in the rest of the unit cells being their translated copies. This fact is a 
consequence of translation symmetry and can be stated mathematically as 

|a(R i ))=T(R i )|a(0)), (1) 

where |a(0)) represents a Wannier orbital localized in the reference unit cell assumed to be 
located at the origin while |a(Rj)) is the corresponding orbital of the i-th unit cell located 
at position Rj, where Rj is a lattice vector. The corresponding translation is induced by 
the operator T(R,). Clearly, in a Wannier-function-based description, the entire crystal can 
be described by the orbital set {\a(Rj)); a = l,n c ;j = l,N}, where n c is the number of 
Wannier functions associated with a unit cell and N is the total number of unit cells in the 
crystal whose limiting value, of course, is infinity. In the following discussion Greek labels 
a, /?,... , will be reserved for labeling the (occupied) Wannier orbitals of a unit cell. 

Having established the fact that one needs to specify the Wannier orbitals of only one of 
the unit cells, we need a dynamic prescription for obtaining them. In paper ii we adopted 
an "embedded-cluster" picture of an infinite solid to achieve that goal. In this scheme we 
envision the infinite crystal as a central cluster (the reference cell) embedded in the field of 
its environment consisting of an infinite number of self-similar unit cells arranged according 
to the crystal lattice. Within this picture, it is straightforward, for a closed-shell system, to 
write down the restricted Hartree-Fock (RHF) equations for the Wannier functions of the 
embedded cluster, under the Born-Oppenheimer approximation: 

(T + U + YpJp A$| 7 (Rfc)><7(R*)|)|a> = e a |a>, (2) 

keAf 7 

where |a) stands for |a(0)), an orbital centered in the reference unit cell, T represents the 
kinetic-energy operator, U represents the interaction of the electrons of the reference cell 
with the nuclei of the whole of the solid, while Jp, Kp are the conventional Coulomb and 
exchange operators defined as 



Ma) = E#(R,)I^I/W>I«> 
Kp\a) = E,(/3(R J )I^I«)|/5(R,)) 



(3) 



Although here we have outlined an intuitive derivation of Eq.(§), it can alsobe obtained 
rigorously by minimizing the total energy of the infinite solid as shown in IlE-3. The first 
three terms on the left hand side of Eq. (0) constitute the conventional canonical Hartree-Fock 
operator while the last term plays the role of a localization potential. It involves a sum over 
projection operators constructed from the orbitals of the unit cells localized in the immediate 
neighborhood M of the central cluster. For an infinitely high value of the shift parameters, 
A*, it has a localizing effect on the orbitals of the reference cell while simultaneously making 
the latter orthogonal to those located in Af. The choice of Af will clearly be dictated by the 
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localization characteristics of the electrons of the system under consideration. It is expected 
to be larger for systems with smaller band gaps. In our calculations we have typically chosen 
M to include up to third nearest-neighbor unit cells of the reference cell. The orthogonality 
of the orbitals contained in unit cells beyond J\f will also be dependent upon the band gap 
of the system and should be automatic once the region M has been chosen to be sufficiently 
large. Numerical values in the range 10 3 — 10 5 atomic units for the shift parameters were 
found to be suitable. 

We have adopted a linear-combination of atomic orbital (LCAO) formalism utilizing 
lobe-type Gaussian basis functions^ to solve Eq.(^). Terms U, J and K appearing in Eq.(^) 
involve infinite lattice sums and deserve special consideration. Evaluation of these terms 
along with other computational aspects are discussed in detail in paper It is clear 

that the orthogonalization of the orbitals of the reference cell to those of the region Af will 
introduce oscillations in these orbitals which are also referred to as the orthogonalization 
tails. In order to describe these orthogonalization tails, one needs to express the orbitals 
of the reference cell as linear combination of basis functions located both in the reference 
cell as well as in N ' . This increases the dimension of the Fock matrix to be diagonalized 
as compared to a canonical Hartree-Fock program as implemented, e.g., in crystal!. 
This, however, does not affect the evaluation of one- and two-electron integrals (and their 
number) if a careful use of translational symmetry is made as discussed in 110. After the 
diagonalization of the Fock matrix, one is confronted with the task of choosing the occupied 
orbitals from the spectrum of the eigenvalues. For the ionic systems considered here, the 
aufbau principle was used for this purpose; in practice, one starts with a suitable guess 
for the Wannier functions, and iteratively solves Eq.(|2]) until the energy per unit cell has 
converged. 



III. CALCULATIONS AND RESULTS 

In this section we present the results of the calculations performed on crystalline Li 2 
and Na20. These compounds have been investigated extensively at the Hartree-Fock level in 
a number of papers by the Torino group using their CRYSTAL programE^-a. Therefore, we 
also intend to compare our results with the most recent results of that groupB In addition, 
we discuss the characteristics of the O 2- Wannier orbitals obtained from our calculations 
and compare them to atomic orbitals of free oxygen (and its singly negative ion). 



A. Geometry, basis set and computational parameters 

All the calculations reported below have been carried out by assuming the observed anti- 
fluorite structure with the space group Fm3m. The reference unit cell was assumed to be 
the primitive cell with oxygen anion placed at the origin and the two cations at the positions 
(a/4, a/4, a/4) and (—a/4, —a/4, —a/4), where a is the lattice constant. 

The Torino group studied these compounds using highly-optimized extended basis sets 
and for the sake of comparison, we have used the most recent basis set reported by thenfi It 
employs a (7slp)/[2slp] set for lithium, a (15s7p)/[4s3p] set for sodium and a (14s6p)/[4s3p] 
set for oxygenEJ. As indicated in the previous section, we have not used real Cartesian-type 
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Gaussian basis functions but rather their approximate counterparts obtained by forming 
linear combinations of lobe-type (Is) Gaussian functions0@. For the case of Li 2 the 
number of (contracted) basis functions per unit cell is 23 while for the case of Na 2 it is 
39. However, in order to satisfy the orthogonality requirements associated with the Wannier 
functions, we also have to supply the basis functions located on the neighboring cells to 
describe the Wannier orbitals of the central cluster. As mentioned in the previous section, 
this neighborhood M consists of up to third-nearest neighbors of the central cluster. Since 
for the fee geometry, there are 42 unit cells in M so defined, one has to use the basis functions 
of 43 unit cells including the reference cell. Thus the number of basis functions associated 
with the Wannier orbitals of Li 2 and Na 2 explodes to 989 and 1677 respectively. This 
proliferation in the number of basis functions does not affect the integral evaluation time in 
any way0, compared to the CRYSTAL program!, as long as the use of point-group symmetry 
is the same. However, it does increase the dimension of the Fock matrix to be diagonalized 
to the corresponding numbers, thus increasing the time needed to perform the Hartree-Fock 
iterations. The latter problem is not too critical, though, since contracting basis functions in 
Af to about single-zeta quality is expected to remedy this issue without adversely affecting 
the accuracy of the calculation. 

The CRYSTAL program uses several computational parameters which determine the 
accuracy of the Coulomb and the exchange series. The parameter related to the Coulomb 
series is called ITOL1 and those related to the exchange series are called ITOL3, ITOL4 
and ITOL5& The values of 7,7,7 and 15 for these parameters for the CRYSTAL-program 
based calculations are generally believed to lead to well-converged results!. These values of 
the parameters ensure an absolute accuracy of ~ 1.0 x 10~ 7 atomic units (a.u.) in the Fock 
matrix elements leading to an expected accuracy of ~ 1 milliHartree per atom in the total 
energyi. Therefore, to make the comparison with CRYSTAL results transparent, in our 
calculations also we treated the Coulomb and the exchange series in an entirely equivalent 
way. 



B. Lattice Constant, Bulk Modulus and the Cohesive Energy 

In order to optimize the lattice constant and to obtain the bulk modulus we first com- 
puted total energies per unit cell of both the compounds for different values of lattice con- 
stants which are closely spaced around the equilibrium value. As also reported in papers 
I and II§0, the total energies (per unit cell) obtained by us at different lattice constants 
agreed with the corresponding values obtained by the CRYSTAL program to within frac- 
tions of a milliHartree. We then fitted these data points to polynomials. The resulting 
bulk properties were found to be stable to within a fraction of a percent with respect to 
an increase in the degree of the polynomial considered, and are listed in table |. In the 
same table, lattice constants, bulk moduli, and cohesive energies for the two compounds 
obtained from our calculations are compared with the results of the Torino groupB and also 
with experiments whenever possible. Cohesive energies were obtained by subtracting the 
corresponding atomic Hartree-Fock energies from the equilibrium total energy per unit cell. 
To make the comparison with the results of the Torino group meaningful, we used the same 
atomic HF energies of -7.4313 a.u. (Li atom), -161.8513 a.u. (Na atom) and -74.8012 a.u. 
(O atom) as used in their calculations^. 
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With the same basis sets, our approach and that of the Torino grou p0 should yield 
identical results. A quick glance at table [j] reveals that the agreement for lattice constants 
and cohesive energies is excellent. The maximum disagreement of 1 milliHartree in the 
cohesive energies is well within the expected numerical accuracies of crystal! and our 
program. The maximum deviation in the bulk moduli is less than 3 GPa and are to be 
expected for the present disagreements in total energies, since second derivatives are much 
more sensitive to numerical inaccuracies. We suspect that part of the reason behind these 
small differences in the results of the two programs could also be due to our use of lobe 
functions to approximate the Cartesian-type Gaussian basis functions used in the CRYSTAL 
program!. 

The experimental value of the lattice constant of Li20 is based on inelastic neutron 
scattering experiment sEil performed in a temperature range 293-1603 K. The value of 4.573 
A is the zero-temperature value obtained by extrapolating the a versus T curve@. The 
result of the Torino erou is in essentially exact agreement with this value while our value 
is 0.003 A shorter. For Na20 the experimental value of the lattice constant was determined 
to be 5.55 A at room temperature by Zintl et al.0 in an old experiment based upon powder 
diffraction pattern data. Because of the lack of availability of any other measurements, 
Dovesi et al.0 extrapolated this to a zero-temperature value of 5.49 A using arguments 
based upon trends observed in the temperature dependence of the lattice constant of NaF. 
Our calculated value of 5.481 A as well as the value of 5.484 A of the Torino group are 
both in good agreement with the above-mentioned zero-temperature value. 

The extrapolated zero-temperature value of the bulk modulus of Li20 obtained from the 
finite temperature data of Hull et al.l^l is 89 GPa. Our value of 94.6 GPa is approximately 
6% larger than this experimental value. We believe that most of this disagreement is due to 
missing correlation effects. In the case of Na20, no experimental value is available for the 
bulk modulus, to the best of our knowledge. However, reasonably close agreement with the 
results of the Torino grou gives us confidence as to the correctness of our result. 

When we compare the Hartree-Fock cohesive energies reported here to the experimental 
values0, we note that for Li 2 we recover « 67% of the total contribution while for Na20 
this fraction is about 57%. As also argued by Dovesi et al.0, most of the missing cohesive 
energy at the Hartree-Fock level is due to the absence of correlation effects. This belief is 
also substantiated by various finite-cluster-based calculations performed in our group where, 
by including correlation effects, typically about 95% of the experimental cohesive energy was 
recoveredta. Therefore, the extension of the scheme of local correlation-energy increments 
as described in refs.@, to the case of infinite crystals, appears to be worthwhile. This task, 
however, is nontrivial and presently we are in the process of implementing it. 



C. Wannier Functions 

Now we turn to the discussion of the Hartree-Fock Wannier functions that we obtain on 
solving Eq.(Q). Since the qualitative behavior of the Wannier functions of both compounds is 
expected to be the same, we will discuss only the orbitals of Li20. The central cluster chosen 
in this work consists of one primitive cell of the lattice. Therefore, for Li20 our Hartree- 
Fock equations involve fourteen electrons leading to seven Wannier functions. Out of these 
seven functions, three correspond to low-lying ls-type core orbitals of the three atoms in 
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the unit cell, while the remaining four correspond to the 2s- and 2p-type functions centered 
on oxygen. Since these four high-lying valence Wannier functions are responsible for most 
of the chemical properties of the compounds considered, we will restrict our discussion to 
them. 

The nearest-neighbor environment of the O ion in the crystal consists of Li + ions lo- 
cated at the eight corners of a cube with coordinates (±a/4, ±a/4, ±a/4), while the center of 
the cube at position (0,0,0) is occupied by O ion itself. This type of environment results in 
three 2p-type Wannier functions of the crystal, each one of which is a mixture of the 2p x , 2p y 
and 2p 2 type basis functions localized on oxygen. In Figs. [I] and £| respectively, we present 
the 2s-type and one of the 2p-type Wannier functions of Li 2 crystal, localized on the O 
ion. The 2p-type Wannier function presented here is isotropically oriented along the body 
diagonals. Both the orbitals are plotted along the crystal [111] direction. Since the crystal 
structure of Li 2 is invariant under the parity operation, we expect the corresponding 2s- 
and 2p-type Wannier functions to be, respectively, symmetric and antisymmetric, under the 
operation of parity. This is precisely what we find when we examine those figures. In addi- 
tion, we also note that both the orbitals have nodes near the positions (±a/4, ±a/4, ±a/4), 
which correspond to the locations of the two Li atoms of the unit cell. These nodes are a 
consequence of the orthogonalization of the Wannier functions to the Li ls-type core Wan- 
nier functions. The localized nature of both the Wannier functions is evident by their rapid 
decay as one moves away from the oxygen site. This is a rather pictorial confirmation also 
of the ionic character of the crystal. The noteworthy point is that apart from the presence 
of the orthogonality nodes which are indicative of the presence of the environment, the qual- 
itative nature of Wannier functions is very similar to those of any molecular orbitals, thus 
bringing them close to the intuition of a chemist. 

The simplest qualitative picture of cohesion in an ionic solid is that it is accompanied by 
charge transfer from one set of atoms to another, leading to cations and anions at different 
sites. These oppositely charged ions, in turn, are held together by the electrostatic attraction 
between them. For the present case of Li20, this would imply that the crystal consists 
of Li + and O ions in closed-shell configurations. Such a picture can be verified, e.g. 
by performing a Mulliken population analysis, also in a Bloch-orbital-based approachEil. 
However, a Wannier-function-based approach provides a more direct and unambiguous view 
of the cohesive process. To this end we first compare the 2s orbital of the free oxygen atom 
(which is also plotted in Fig. [I]) with the 2s- type Wannier function of Li20. For free oxygen 
atom a basis set also reported by the Torino group0 was used, which they obtained by 
reoptimizing the outermost exponents of the crystal basis set and by adding one more sp- 
type diffuse exponent (0.0867) to it. As is clear from the figure, the spatial behavior of the 2s 
orbital changes only in a rather subtle manner in that the corresponding Wannier function 
develop nodes near the positions (±a/4, ±a/4, ±a/4) due to their orthogonalization to the 
ls-type Wannier functions of the Li + ions of the unit cell located there. Apart from these 
orthogonality nodes, the two orbitals are remarkably similar. This is consistent with our 
intuitive picture in that being a lower valence orbital, one should not expect the 2s orbital to 
be affected appreciably by the crystal field. Therefore, we expect the effects of cohesion to be 
most evident in the case of the 2p-type Wannier functions. In the simple ionic picture, the 
O ion of the solid has a filled 2p shell. The presence of two extra electrons on oxygen will 
result in the 2p electrons experiencing larger on-site repulsion leading to the charge density 
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associated with 2p-type Wannier functions becoming more diffuse compared to the isolated 
atom. Using the same intuitive reasoning we expect the 2p charge density in the solid to 
be more diffuse even compared to the free O . These points are illustrated in Fig. where 
we plot the total 2p charge density associated with the 2p-type Wannier functions of the 
Li20 lattice along the [111] crystal direction, and compare it with the corresponding charge 
densities of free O atom and free 0~ ion. The orbitals of free 0~ ion were obtained from a 
basis set which was obtained by augmenting the free atom basis set discussed earlier by one 
more diffuse sp-type exponent of 0.0345. The total normalized charge densities were obtained 
by proportionately adding the contributions of the individual charge densities of the three 2p- 
type orbitals of the systems under consideration. A derealization trend is certainly obvious 
in the region close to the nuclei from the decreasing heights of the peaks as one compares 
the 2p-charge densities in the order of free O atom, free 0~ ion and embedded O ion. To 
quantify the degree of localization of various orbitals we computed the directionally averaged 
expection values of r 2 = x 2 + y 2 + z 2 operator for the 2p-type orbitals of free O, 0~ and the 
corresponding Wannier functions of the Li 2 crystal. The directionally averaged values of 
a given system were obtained by averaging over the expectation values of the three 2p-type 
orbitals according to < r 2 > 2p = (J2i=i < r2 >2 P J/3.The obtained values for < r 2 > 2p (in 
Bohr 2 ) were 1.974 (free O atom), 3.038 (free 0~ ion) and, 3.008 (Li 2 crystal). As expected, 
we find that the value of < r 2 > 2p is smallest for the case of free O atom. However, when 
we compare the corresponding values of free O - with those of the embedded O of the 
Li20 crystal, we arrive at somewhat counter-intuitive result that the size of the 2p-type 
orbital of the embedded O is a little smaller than the corresponding orbital of free O . 
To account for finite-basis-set effects as well as the nonuniqueness of Wannier functions, we 
can perhaps say that the 2p-type orbitals of free 0~ and embedded O are approximately 
equal. One possible objection against this finding is that our basis set for the solid-state 
calculations, with the most diffuse sp-type exponent of only 0.150, is not diffuse enough 
to describe a doubly charged anion like O . However, we do not suspect this to be the 
reason behind our finding because, as discussed earlier, our Wannier functions contain basis 
functions not only of the atom that they are centered upon, but also the basis functions of 
atoms contained in up to third-nearest neighbour unit cells of the reference cell. Thus these 
moderately diffuse exponents located on the neighboring atoms should be able to describe 
the long-range behavior of the Wannier functions of the reference atoms. Moreover, since the 
solid state basis set has correctly predicted the lattice constant of Li 2 0, we expect it to be a 
reasonably complete set. Our explanation of the finding, that the size of embedded O ion 
of the Li 2 crystal is essentially the same as that of the singly charged free 0~, is based upon 
a competition between the on-site repulsion and the crystal-field effects (Pauli repulsion of 
surrounding atoms). In our opinion, in the crystalline phase, the on-site repulsion which 
certainly has a delocalizing influence, is overpowered by the crystal field leading to a stable 
O ion whose size is comparable to that of the free O - ion. 



IV. CONCLUSIONS 

In this paper we extended the application of a Wannier-function-based Hartree-Fock 
approach developed earlieri'0 to alkali oxides, Na 2 and Li 2 0. State-of-the-art basis sets 
containing rather diffuse s- and p-type functions were used and quantities such as the equi- 
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librium lattice constant, bulk modulus and cohesive energy per unit cell were computed. We 
obtained excellent agreement with the results of a previous Hartree-Fock study carried out 
by the Torino groupB using their Bloch-orbital-based CRYSTAL program!. In addition, a 
detailed pictorial view of the cohesive process is provided by the Wannier functions of our 
work which comes close to the intuitive understanding of a chemist. 

As far as the agreement with experimental data is concerned, it is excellent for the 
lattice constant. For the bulk modulus agreement is acceptable. As expected, the most 
severe deviation is evident in case of the cohesive energies of the two compounds where 
disagreement between Hartree-Fock theory and experiment is in the range of 30% — 40%. 
This points to the importance of correlation effects, and we consider it as an advantage of 
our Wannier-function-based Hartree-Fock approach that it provides a natural starting-point 
for improvement by means of local correlation methodsEl. Efforts along this direction are 
presently underway in our group, and the results will be presented in a future publication. 
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FIGURES 



FIG. 1. Free atom 2s orbital of oxygen and 2s-type crystal Wannier function of Li20 localized 
on O plotted along the [111] direction. 

FIG. 2. 2p-type crystal Wannier function of Li20 localized on O plotted along the [111] 
direction. 

FIG. 3. Charge density associated with the 2p-type Wannier functions of Li20 localized on 
O , compared to the 2p charge densities of the free oxygen atom and free O . All the charge 
densities are plotted along the [111] direction and are normalized to unity. 
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TABLES 



TABLE I. Equilibrium lattice constants (in A), bulk moduli (in GPa), and cohesive energies 
(in atomic units) for Na20 and Li20, obtained using our approach and reported by the Torino 
erou Relevant experimental data are also given for comparison. 



Quantity 


Method 


Li 2 


System 


Na 2 




This Work 


4.0 1 u 






Lattice Constant 


Torino 


4.573 




5.484 




Exp 


4.573 a 




5.55 b 
(5.49) 




This Work 


94.6 




61.1 


Bulk Modulus 


Torino 


92.6 




58.4 




Exp 


89 a 








This Work 


-0.3008 




-0.1883 


Cohesive Energy 


Torino 


-0.3005 




-0.1893 




Exp 


-0.4491 c 




-0.3321 c 



a RefM 



Room temperature value from ref.Ha The number in parentheses is the extrapolated zero- 
temperature value suggested by the Torino group0. For a discussion see text. 
c RefM 
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